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ABSTRACT 

Half a year after its outburst in September 2002, nova V4 743 Sgr evolved into the brightest 
supersoft X-ray source in the sky with a flux maximum around 30 A. We calculated grids of synthetic 
energy distributions (SEDs) based on NLTE model atmospheres for the analysis of the hottest 
white dwarfs and present the result of fits to Chandra and XMM-Newton grating X-ray spectra of 
V4743 Sgr of outstanding quality, exhibiting prominent resonance lines of C v, C vi, N vi, N vii, and 
O VII in absorption. The nova reached its highest effective temperature (Tcff = 740 ± 70 kK) around 
April 2003 and remained at that temperature at least until September 2003. We conclude that the 
white dwarf is massive, « 1.1 — 1.2 M©. The nuclear-burning phase lasted for 2 to 2.5 years after 
the outburst, probably the average duration for a classical nova. The photosphere of V4743 Sgr was 
strongly carbon deficient (« 0.01 times solar) and enriched in nitrogen and oxygen (> 5 times solar). 
Especially the very low C /N ratio indicates that the material at the white dwarf's surface underwent 
thermonuclear burning. Thus, this nova retained some of the accreted material and did not eject all 
of it in outburst. From March to September 2003, the nitrogen abundance is strongly decreasing, 
probably new material is already been accreted at this stage. 
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1. INTRODUCTION 

Nova V4743 Sgr (Nova Sgr 2002 c) was discovered in 
outburst in September 2002 and reached V ~ 5 mag on 

2002 September 20 (Haseda et al. 2002). It was a very 
fast nova, with a steep decline in the optical light curve 
and large ejection velocities. The time to decay by 3 mag 
in the visual (ta), was 15 days and the FWHM of the 
Ha line reached 2400kmsec-i (jKato et al.l 12002^. Es- 
timates of the dista nce from infrared observat ions vary 
from 1 200±300pc ( Nielboc k fc Schmidtobreickll2003D to 
« 6 300pc (Starrfield & Lyke priv. comm.). 

In December 2002, the nova was observed for the first 
time with Chandra ACIS-S. At that time it was a very 
soft and moderately luminous X-ray source, with a count 
rate of about 0.3 cts sec~^ . There were indications that it 
was not at the peak of X-ray luminosity yet, so a Chan- 
dra HRC-fLETG observation was done only later, on 

2003 March 20. The count rate was astonishingly high, 
40 cts sec^^ during 3.6 hours, then a slow decay, last- 
ing for an hour and a half, was followed by another 1.5 
hours of very low luminosity with a measurement of onl y 
« 0.02 cts sec-i (|Ness et ai.ll2"00l IStarrfield et al.ll200l . 

|rauch@astro.uni-tuebingen.de, orio@astro.wisc.edu| 



This decline was not due to an eclipse, since the orbital 
period of the syste m is 24 278 ± 259 sec (or 6.74 hours, 
IWagner et al.|[2003l ) and no eclipse was observed in a fol- 
lowing 10 h XMM-Newton observation two weeks later. 
An obscuration of a supersoft X-ray source in a nova 
was observed on ce before with BeppoSAX, in V382Vel 
(|Orio et al.ll2002[ ). The only reason for this sudden nearly 
shut down of the source may have been the ejection of 
a new shell of material that w as optically thick to su- 
persoft X-rays from the surface (|Orioll2008D . During the 
Chandra observation in March of 2003 an oscillation with 
a period of 1 324 sec was detected, with fluctuations of 
20% from the me an count rate (see Starrfield et al. 200j; 
INess et al.|[200l . The most likely cause of this oscilla- 
tion is a non-radial pulsation of the WD (Starrfield et al.l 
I2OOI . 

A second observation was proposed by us to the XMM- 
Newton Project Scientist as a target during the Discre- 
tionary Time. The nova w as observed for 1 h with this 
satehite on 2003 April 4 (jOrio et al.|[200l . Three X- 
ray telescopes with five X-ray detectors were all used: 
the E uropean Photon Imaging C amera pn (IStriider et al] 
l200ll ). two EPIC MOS cameras ([Turner et al.ll2001[ ). and 
two Reflection Grating Spectrometers (RGS-1 and RGS- 
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2, iHerder et al.|[200ll ). The observation lasted a little 
over 36 000 sec. The data obtained by the two EPIC 
MOS cameras suffered from very severe pile up, but the 
operation of the pn camera was switched to timing mode 
and yielded a measured average, backgroun d corrected 
EPIC-pn count rate of 1 348.0 ± 3 cts sec"^ (lOrio et all 
I2003D . with variations by 40% (jLeibowitz et all I2006D . 
The RGS count rates were about 57 cts sec~^, and the 
unabsorbed flux was 1.5 x 10~^ erg cm~^ sec~^, consis- 
tent with the flux measured 16 days earlier with Chandra 
(jLeibowitz et al.|[200 6). Most of the variability is well 
represented as a combination of oscilla tions at a set of 
discrete frequencies lower than 1.7 mHz (jLeibowitz et al.] 
|2006( ). At least five frequencies preserve their coherence 
over the 16 days time interval between the two observa- 
tions. The 1324 sec period detected with Chandra could 
be split into two periods, of 1310.1 sec and 1371.6sec, 
both of which are consistent i n the Chandra data b ut it 
may not have been resolved (jLeibowitz et al.l 120061 ). In 
that article we suggested that a period in the power spec- 
trum of both light curves at the frequency of ~0.75mIIz 
(corresponding to 1371.6 sec) is the spin period of the 
white dwarf in the system, and that other observed fre- 
quencies are signatures of non-radial white dwarf pulsa- 
tions. 

The X-ray evolution was followed with 10 000 sec long 
exposures with the Chandra LETG in July and Septem- 
ber 2003 and in February 2004, and with a 30 000 sec 
observation with XMM-Newton in September 2004. By 
this time, the nova X-ray luminosity had decayed by 5 
orders of magnitude, but the period of w 1370 sec was 
still detectable. However, the S/N of the RGS spectra at 
this epoch was very poor. 

We have performed a NLTE spectral analysis of the 
grating spectra obtained between 2003 March and 2004 
February, with NLTE model-atmosphere techniques. 
The TMAP model-atmosphere code and the construc- 
tion of the model atoms that are used are described in 
Sect. 12. l1 An attempt to determine photospheric parame- 
ters with a XSPEC fit procedure is presented in Sect. 13. II 
This is followed by a detailed investigation of the ob- 
served, flux-calibrated spectrum of V4743Sgr f Sect. 13. 2) ). 
In Sect.m we describe how our models fit the Chandra 
and XMM-Newton observations to interpret the post- 
nova evolution in the 18 months after the outburst of 
V4743Sgr. 

2. THE USE OF MODEL ATMOSPHERES 

A previous analysis of Chandra LE TG-S obser- 
vations of V4743Sgr (jPetz et al.l I2005D . done with 
the PHOENIX NLTE models with solar abundances, 
reached Teg = 580 kK. The fit of the model to the data 
was still poor, and this is not surprising because in the 
phase directly after the outburst, i.e. when the H burn- 
ing is still on-going (for years), the surface composition 
of the WD is poorly known, but it is very unlikely to 
be solar. With the PHOENIX code, the nova atmo- 
sphere is approximated as a n expanding, but stationary 
in time structure. Recently, Ivan Rossum fc Ne ss (20lJ) 
presented a new version of PHOENIX that considers 
mass-loss as well as velocity fields. In the future, such 
codes will become a powerful tool for the analysis and 
understanding of novae and other supersoft sources. 

However, in order to make progress, we decided to 



neglect the velocity field, and to use our static NLTE 
models to investigate basic parameters like Tog and sur- 
face composition. This is not fully justified because es- 
pecially in the first spectra the lines were significantly 
blue-shifted, but our aim is at least a qualitative model- 
ing of the XMM-RCS and Chandra-LETC observations 
done between 2003 March and 2004 February. We started 
with the highest S/N spectrum, the one obtained with 
the XMM-Newton RGS gratings on April 4 2003, that 
we used as a template to adjust the abundances. We 
then proceeded to fit also the Chandra spectra with the 
same model, checking whether the abundances we ob- 
tained were suitable, and adjusting T^g for each epoch. 
In this way, we obtained the evolution of Teff, and the 
duration of the constant bolometric luminosity phase. 




wavelength / A 

Fig. 1. — Comparison of synthetic spectra (Tcff = 700kK, 
logg = 9), from models with different elemental composition (see 
labels). The abundances are given in Table[3] (model A). For clar- 
ity, the HHe and the H-Fe spectra are shifted in log Fx by -\-2 and 
—4, respectively. 



2.1. Model atmospheres and atomic data 

For a reliable analysis of X-ray observations of the 
hottest white dwarfs, detailed NLTE model-atmospheres 
that consider opacities of all elements from hydro- 
gen up to the iron-group elements (cf. iRauchI l2003f) 
are required. Thus, for our analysis, we employed 
the plane-parallel, static models calculated with the 
Tiibingen NL TE Model- Atmosphere Package {TMAI^, 
I Werner et al.l l2003). The constructio n of model atoms 
which are used within TMAP follows iRauch fc DeetienI 
(j2003l ). Some details for these extremely hot model at- 
mospheres in the MK-Toff range are summarized briefly 
in the following. 

Since the surface composition is unknown, we started 
with the calculation of exploratory H-|-He-|-C-|-N-f O 
models and later included Ne, Mg, Si, S, as well as 
the iron-group elements. The iron-group elements (Sc 
- Ni) and Ca are trea ted with a statistical method 
(jRauch fc DeetieiJ l2003f) and are represented by one 
generic model atom. 

Beside an element selection, also the construction of 
the model atoms has to be performed with care. An 
unrealistic upper cut-off in the number of considered 

1 http://astro.uni-tuebingen.de/~rauch/TMAP/TMAP.html 



iron-group ionization stages causes an artificial over- 
population of the highest ionization stage, and thus, af- 
fects its lines and the flux level. In Figure[Tl the impact of 
iron-group opacities on the astrophysical flux is demon- 
strated. 
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Fig. 2. — Ionization fractions of the generic iron-group (IG) model 
atom in a T^ff = 700kK, log 3 = 9 model. IG XVII is dominant in 
the line-forming region (—5.6 < logm < 0.5). 



For all values of Tcs, the necessary ionization stages 
of all atoms are determined in advance by test calcula- 
tions. The selection criterion is, that at least the ion 
(e.g. in the case of species X), X"+ that is dominant in 
the line-forming region has to be included together with 
the neighboring two, i.e. X("-i)+ and X("+i)+. Thus, 
the model atoms contain in general three to five ioniza- 
tion stages. E.g. in the case of our generic iron-group 
model atom, IG xvii is dominant at Tcff = 700kK and 
log 3 = 9 [cm sec~^] (Figure[2]). We therefore selected the 
ionization stages xiv - xviii. Statistics of the model 
atoms are shown in Table[T] In total, 228 atomic levels 
are treated in NLTE, 360 additional levels in LTE, and 
349 individual line transitions are considered. 

3. SPECTRAL ANALYSIS 

A first grid of models is co mposed of H-|-He-|-C -|-N+0 
with solar abundance ratios ([Asplund et al.|[2009() within 
Tcff = 0.45 — 1.05 MK and a fixed surface gravity of 
log.g = 9 (Figure[3|). We note that all synthetic energy 
distributions (SEDs) in our model grids described here 
are available at in Virtual Observatory ( VO) compliant 
form from the VO service TheoSSM provided by the 
German Astrophysical Virtual Observatory {GAV(^) as 
well as atablefS for the use with XSPE(M 

3.1. Preliminary analysis 

We started with the XMM-Newton RGS spectra of 
2003 April (TableE]). The data were reduced with 
the ESA XMM-Newton Science Analysis SystenQ (SAS) 
software, version 5.3.3, using the latest calibration files 
available. The RGS dispersion gratings cover the 5 — 35 A 
wavelength range (2.48 — 0.35 keV), although we obtained 

^ http:/ /vo. ari.uni-heidelberg.de/ssatr-0.01/TrSpectra.jsp? 
^ http://www.g-vo.org 

* http://astro.uni-tuebingen.de/~rauch/TMAF/TMAF.html 
^ http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec 
® http://xmm.esa.int/sas/ 
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TABLE 1 

Statistics of the model atoms used in our calculations of 
THE NLTE model atmospheres with Tcff = 0.7 MK. The 

NOTATION IS: N = LEVELS TREATED IN NLTE, L = LTE LEVELS, 
R = RADIATIVE BOUND-BOUND TRANSITIONS. FOR THE GENERIC 
IRON-GROUP (IG) MODEL ATOM (WHICH REPRESENTS Ca, SC, Ti, 

V, Cr, Mn, Fe, Co, and Ni) the number of individual lines 

WHICH ARE combined TO SO-CALLED SUPERLINES IS GIVEN IN 
BRACKETS. 
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FiG. 3.— Comparison of NLTE model fluxes (H+He-|-C-|-N-|-0, 
solar abundances, top panel) at different T^g. The wavelength 
interval corresponds to 0.95 — 0.26 keV. The positions of C, N, and 
O lines (at top) and bound-free absorption thresholds are marked. 
In the bottom panel, the abundance ratios are taken from nova 
models of Prialnik (priv. comm.) except for the carbon abundance 
which is reduced by a factor of 100. 

useful signal only in the 5 — 26 A range (2.48 — 0.48 keV). 
The RGS data are piled-up, but in a dispersion instru- 
ment, piled- up events at a discrete wavelength increase 
in pulse height amplitude by an integer multiple of the 
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intrinsic energy. Furthermore, since the softness of the 
source precludes any intrinsic photons from higher spec- 
tral orders, we can confidently identify events that occur 
within the higher order spectral masks for an on-axis 
point source as piled-up first order photons. This is ver- 
ified by line matching of the piled-up events using the 
first order response matrix. Source events dominate over 
background and scattered source light in the first three 
orders. Consequently we added events within the sec- 
ond and third order extraction masks to the first order 
events, thus reclaiming the piled-up events and increas- 
ing the signal-to-noise of the spectra. 

In a first step, we compared our SEDs with the XMM- 
Newton observation of V4743Sgr via XSPEC. We let 
the neutral hydrogen column density TVhi [cm~^] and 
Teff vary as free parameters in XSPEC. XSPEC deter- 
mines Toff = 709 kK and logiVni = 20.58. In Figureg] 
(top panel), we show this XSPEC fit. The reader may 
miss the typical XSPEC residuals at the bottom of the 
plots. Since our models do not include all the elements 
with all the lines that may be exhibited in the observa- 
tion, the residuals are not that helpful like in comparisons 
of simple models where continuum slope and a handful 
of strategic lines has to be reproduced. Thus, we decided 
to omit the residual plots generally. Instead, the reader 
should judge the strengths of prominent lines and ab- 
sorption edges with respect to the (estimated) local flux 
continuum. 

As expected, a solar composition does not well repre- 
sent the white dwarf spectrum in V4743 Sgr. The flux 
at wavelengths smaller than 26 A is far too low. In 
test calculations, we decreased the carbon abundance 
to [C] = — 1.0 ([X] denotes log [mass fraction / solar 
mass fraction] of species X) and to [C] — —2.0. We 
achieved a much better fit at [C] = — 2.0, Teff = 610kK, 
and log A^Hi = 20.54 for A < 26A (Figuregl). 

We wondered whether the small wavelength interval 
(18 - 38 A) of our XMM-Newton RGS spectra leads in 
the XSPEC fit proc edure to a smaller iVni than found 
bv lPetz et all (l2005l . logiVni = 21.60, Tefi = 580kK) ob- 
tained fitting the a Chandra LETG spectrum (18 — 58 A). 
In Sect.H] we show however, that we fitted the Chan- 
dra LETG spectrum of M arch 19, 20 3 and used the 
same wavelength range as iPetz et all ()2005[ ). deriving 
Teff = 601kK and logiVni = 20.95 for the [C] = -2.0 
model. A decreasing Nhi from March to April is con- 
sistent with intrinsic absorption of the ejecta clearing 
up. 

3.2. Spectral analysis of the flux- calibrated spectrum 

Before we start now with a detailed analysis, we men- 
tion that the calibration of the instruments that were 
used to obtain our data is still an issue. Different calibra- 
tion products (response matrices and effective areas) may 
play an important role on the results obtained from fits 
to X-ray spectra. Cross calibration be tween instruments 
has improved over the years (see, e.g. IBeuermann et ahl 
[2006, 2008) but there is yet room for further improve- 
ment. 

While the wavelength scale of X-ray gratings instru- 
ments is generally very well known, the effective areas of 
these instruments are less well calibrated. For example, 
a comparison between the continuum calibrations of the 



energy / keV 
0.35 0.41 0.50 0.62 




30 

wavelength / A 



Fig. 4. — Comparison of the XMM-Newton RGS-1 observa- 
tion of V4743 Sgr (TableEJ with NLTE model-atmosphere fluxes 
(H+He+C+N+O models with solar abundances, but different car- 
bon abundances: top: [C] = 0.0, middle: [C] = — 1.0, bottom: 
[C] = -2.0). Toff and A^h are determined with XSPEC. The SEDs 
are shifted by —1800 km/sec in order to match the observed posi- 
tion of the N VII resonance line. The positions of resonance lines 
of C VI, N VI, N VII, O VII, and O VIII are marked. 



gr atings and CCD i nstrum ents on board XMM-Newton 
bv lStuhlinger et all (|2008| ) shows deviations between the 
RGS and the EPIC-pn of 3% above 0.5 keV and about 
10 % below that energy. While not negligible, these devi- 
ations are still smaller than the deviations between data 
and model caused by simplifications in the atmosphere 
modeling, and can therefore be ignored. 

The preliminary analysis presented in the last section 
shows clearly that the carbon abundance in V4 743 Sgr 
is strongly subsolar. The determination of T^g and a 
more precise abundance determination is hampered by 
the complex XSPEC fitting procedure and the necessity 
to calculate extended grids of models with different Tcff 
if the abundances are changed. The latter results in un- 
reasonable computational times. 

A more straightforward approach is to use the flux- 
calibrated spectra instead of count-rate spectra. How- 
ever, these spectra contain uncertainties in the flux. We 
tried to calibrate the RGS spectra of V4743 Sgr in two 
different ways. One way is by making use of the provided 
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TABLE 2 

X-RAY OBSERVATIONS OF V4743 SGR USED IN THE ANALYSIS. 



satellite 


instrument 


Obsid 


Obs date (start) 


exp. time 

sec 


Chandra 


HRC+LETG 


3775 


2003-03-19 09:30:12 


27058 


XMM-Newton 


RGS-1 


0127720501 


2003-04-04 22:24:12 


35215 


XMM-Newton 


RGS-2 


0127720501 


2003-04-04 22:24:12 


35214 


Chandra 


HRC+LETG 


3776 


2003-07-18 21:35:29 


13746 


Chandra 


HRC+LETG 


4435 


2003-09-25 23:51:50 


14378 


Chandra 


HRC+LETG 


5292 


2004-02-28 01:07:42 


13530 



energy / keV 
0.35 0.41 0.50 0.62 
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Fig. 5. — Comparison of the RGS fluxed spectrum of V4743 Sgr 
obtained from the XMM-Newton pipeline (thick line) with the ab- 
solute flux calculated by XSPEC (thin line, see text for details). 
The XSPEC spectrum is convolved with a Gaussian of 0.10 A 
(FWHM) in order to match the resolution of the binned RGS fluxed 
spectrum. Note that both spectra agree well for A^32 A. The marks 
indicate positions of lines that are included in our models. 



XMM-SAS pipeline products, the so called "fluxed" 
spectra. Such spectra are e.g. shown by i?ji?_CQ, the 
Browsing Interface for RGS Data, with the purpose of 
visualizing the data free from the peculiarities of the in- 
strument. However, RGS fluxed spectra are computed 
in the pipeline by dividing the extracted spectrum by 
the effective area calculated from its corresponding re- 
sponse matrix. This procedure neglects the redistribu- 
tion of monochromatic response into the dispersion chan- 
nels, and therefore these spectra should not be used for 
very detailed analysis. The second version of the flux- 
calibrated spectrum is calculated by XSPEC and is not 
independent from the models which are used within the 
XSPEC fit procedure (Figure[5|). 

We are aware of these uncertainties, but with the aim 
to determine Tcff within an error range of about 10 % 
and element abundances within 0.5 dex, we used the flux- 
calibrated spectra. We calculated synthetic profiles of 
the resonance lines of the two highest ionization stages 
of C, N, and O (FigureE]). These appear to be strongly 
dependent of Toff within the parameter range of our grids. 

The best XSPEC fit of the [C] = -2.0 model yields 
Tefi = 610kK. Since the maximum fiux level is not 
well reproduced by this fit (Figurel!]), we estimate that 
Teff of V4 743 Sgr may be slightly higher. Thus, we 
selected models with Tcff = 600 — 800 kK and model-A 

^ http://xmm.csac.esa.int/BiRD/ 
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Fig. 6.— Synthetic profiles of the O VII, O VIII, N VI, N VII, 
C V, and C VI resonance lines, calculated from H+He+C+N+O 
model atmospheres with T^ff = 0.4 - l.OMK, logg = 9 and solar 
abundances. The nitrogen lines of the cooler models that appear 
in emission are shifted for clarity. 



abundances (Table[3]) for a comparison with the RGS 
fiuxed spectrum (FigureEj) . The flux level below < 30 A 
is strongly dependent on Toff and at first glance, we 
can achieve good agreement with the RGS fiux level 
at Toff = 700 kK and log TVhi = 20.90. Moreover, the 
prominent absorption lines N vi AA 28.79, 24.90, 23.77 A, 
N VII AA 24.78, 20.91, 19.82, 19. 32 A, and O vii A 21.60 A 
are well reproduced by our model. A detailed inspection 
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of the O VIII A 18.96 A resonance line shows that the ob- 
served continuum flux on its high-energy wing is much 
lower than in the synthetic spectrum. A stronger N vii 
bound-free ground-state absorption edge (18.59 A) could 
probably decrease this flux. Consequently, we increased 
the nitrogen abundance in the following models. 



X x> >x>>xxx >> >> >>xx>x>>>>>>> 
wo ocooococou) CO 2za)c/3zu)OzzzOOO 




35 30 25 20 

wavelength / A 



Fig. 7.— Top: SEDs of our model with [C] = -2.0 and 
Tejj = 600 - 800 kK. They are normalized to the Tg^ = 700kK flux 
level at 38 A for this comparison. The positions of lines are marked. 
Bottom; Comparison of the synthetic spectra of our [C] = —2.0, 
Tcff = 700 kK model with the RGS fluxed spectrum of V4743 Sgr. 



In Figure[7] the Tcft-dependency of the strengths of 
prominent bound-free absorption edges (for an identifi- 
cation, see Figure[9|) is shown. Especially, the N vi and 
N VII ground-state edges at 22.46 A and 18.59 A, respec- 
tively, appear very sensitive on Tog. In order to improve 
the fit, we performed some fine-tuning of the abundances 
(cf. TablelU model B). E.g. in the case of sulfur, model A 
(TablelHl) yields a much too strong S xiv line at 32.51 A, 
i.e. the sulfur abundance is too high in the relevant T^e 
range. In the case of nitrogen, we find prominent absorp- 
tion edges as well as absorption lines in the observation 
which are suitable to adjust the abundance. In the case 
of other species, it is difficult to determine their abun- 
dances precisely because they do not exhibit significant 
features in the XMM-Newton observation. However, this 
can be used to determine at least upper abundance limits. 
Figure[8] demonstrates how the SEDs change if the abun- 
dance of one particular element is artificially increased 
by a factor of ten. We note that only nitrogen, oxygen, 
silicon, and the iron-group elements have an influence on 
the flux in the RGS wavelength range if their abundances 
are increased. 

We compared the strengths of the N vi and N vii edges 
in our models with the observation (Figure[9]). With 



^ ?,><>x X >5>xx>x^^ x^>l^> 
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Fig. 8. — Comparisons of the SED of our Tcfi = 750kK model B 
(thin line, Table|3} to SEDs of models where the abundance of one 
element (indicated by labels) is increased by a factor of ten. The 
positions of ground-state absorption edges are marked at the top. 



TABLE 3 

Element abundances in our models given in [X] = log [mass 

FRACTION / SOLAR MASS FRACTIO N] OF SPECIES X. SO LAR 
ABUNDANCES ARE TAKEN FROM AS PLUND ET AL.I (120091) . We 
NOTE THAT IN MODEL B, THE ABUNDANCES FOR Ne, MG, S, AND 
IG (CA - Nl) ARE DETERMINED UPPER LIMITS. An ERROR RANGE 
OF ±0.3 dex IS ESTIMATED FROM DETAILED COMPARISON WITH THE 
OBSERVATIONS. 



element model A model B 



H 


-0.573 


-0.688 


He 


0.497 


0.382 


C 


-2.614 


-1.513 


N 


0.386 


1.803 


O 


0.425 


1.528 


Ne 


-0.333 


-0.474 


Mg 


0.000 


-0.454 


Si 


0.000 


0.167 


s 


0.000 


-1.583 


IG 


0.000 


0.828 



these models (Tabled model B), we achieve the best fit 
at rcff = 740kK and logiVni = 20.85. 

In FigurefTUl we show a comparison of TMAP SEDs 
with the flux-calibrated Chandra spectra from 2003. In 
2004, the measured flux level and thus. Toff are much 
lower and we disregard that observation here. Further- 
more, for reasons of simplicity, we used a Tcff — 720 kK 
model with model-B abundances (Table[3]) and varied 
only the C and N abundances in order to improve the 
flt. The observed strengths of the Nvi and Nvii reso- 
nance lines as well as their ground-state edges are well 
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energy / keV 
0.35 0.41 0.50 0.62 
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Fig. 9. — Dependance of the strengths of the N VI and N VII 
ground-state thresholds (positions marked by the arrows) on 
Tcff. The SEDs are compared with the RGS fiuxed spectrum of 
V4743 Sgr and normahzed to match the fiux at 38 A. The SEDs 
were convolved with a Gaussian of 0.10 A (FWHM). Positions of 
ground-state thresholds are marked at top, of those line transitions 
just below that marks. 



reproduced for all observations. Since the Nvi / Nvii 
ionization balance is very sensitive to Toff (Figure|6]), 
Teff = 720kK appears a good estimate. Our model fits 
show that the N abundance is decreasing from March to 
September by a factor of about ten. The C lines, e.g. at 
about 34 A are well matched in September while they ap- 
pear too weak in March and July. A higher C abundance 
would result in the appearance of a strong C vi ground- 
state absorption edge at 25.3 A that is not visible in the 
observation. 

In a last step, we used the SEDs of our model-A and 
model-B grids (Table|31) for a comparison with the XMM- 
Newton observation with XSPEC (FigureEI)- The dif- 
ferences in the SEDs of the best fitting models A and B 
are small. If we assume the mean Teff , we are able to con- 
strain the Tcff range of V4743 Sgr to T^q = 740 ± 70 kK. 
The deviations of the abundance of both models (Table[3]) 
show the difficulty of a precise abundance determination 
if only a few spectral lines are available for the analysis. 

4. THE SPECTRAL EVOLUTION OF V4743 SGR 
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Fig. 10. — Comparison of flux-calibrated Chandra observations 
of V4743Sgr (Table^ with our SEDs (Tofi = 720kK, logg = 9^. 
All abundances, but for C and N, are from model-grid B (TablefSj. 
For clarity, the July and September observations are shifted in f\ 
by —2 ■ 10"'"'' and —4 ■ 10"'"'' erg/cm^/sec/A, respectively . 
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Fig. 11. — Comparison of the XMM-Newton RGS-1 observation 
of V4743 Sgr (TableO with our best fitting SEDs from model grids 
A and B (Tabled. and logiVn were determined with XSPEC. 



We then proceeded to fit the Chandra LETG spectra 
at three different epochs using the same models. The 
LETG spectra were reduced using the CIAO softwares, 

^ http://cxc.harvard.edu/ciao/ 
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version 4.1.1. 

In the first 18 months after the outburst of V4743 Sgr, 
Chandra LETG spectra were obtained roughly every 
three months (Table[2]). These observations allow an in- 
vestigation on the spectral evolution, e.g. on the change 
of its Tcff. In Figure[T21 we compare the observations 
with SEDs of our grid (models B). No individual fine- 
tuning of the model abundances for the different spectra 
was performed. Although this might improve the fit, we 
estimate that the impact on the determined T^s is small. 

The Chandra observations (energy range 0.22 — 
0.69 keV) show that V4743Sgr was at a maximum of 
TaS (~ 700 kK) for about one year, then T^s appears 
to decrease by about 10 % in the following six months 
(FigurefT2l left panel). 

In order to judge the quality of our spectral analysis 
of the Apr 2004 XMM-Newton observation that covers a 
smaller energy range (0.32 — 0.69 keV), we restricted the 
fit range of XSPEC (Figure[T21 right panel) for the Chan- 
dra observations. With this restriction, the determined 
Teff is in general higher and the interstellar neutral hy- 
drogen density is smaller, but the deviations are within 
our expected error range of about 10 %. 



5. RESULTS AND CONCLUSIONS 

We have calculated NLTE model atmospheres includ- 
ing opacities of the elements H, He, C, N, O, Ne, Mg, 
Si, S, and Ca - Ni and fitted XMM-Newton RGS spectra 
and Chandra LETG spectra of V4743 Sgr. 

The fit to the RGS spectra (Table[2]) based on these 
models (Table[3l) is best at a T^s of about 740 ± 70kK, 
log 5 = 9, and logA^n = 20.7 ± 0.3. Although this fit is 
not perfect and there are uncertainties in the so-called 
RGS fluxed spectrum used for our spectral analysis, the 
overall flux distribution as well as prominent line features 
are well in agreement with the observation fFigurc fTTj) . 
The photospheric abundances have been adjusted and 
we can give at least upper limits for N, O, Si, and Ca 
- Ni (Figure[51). The determined C/N ratio is very low 
(Table|3]) , indicating that the material on the white dwarf 
surface has been processed through thermonuclear burn- 
ing. We expect that freshly accreted material after the 
outburst or dredged-up material from the white dwarf 
interior, would have a very different range of abundances 
and suggest this is a proof that the white dwarf is re- 
taining some accreted material after each outburst. This 
may imply that V4743 Sgr is on a track towards super- 
nova la explosion, although it retains after each outburst 
only just enough material to burn in less than 2.5 years. 

The fit to the Chandra LETG grating spectra shows 



that V4 743 Sgr reached its highest temperature around 
April 2003 and remained at that temperature at least 
until September 2003. The duration of the constant 
bolometric luminosity phase at constant was be- 
tween 2 and 2.5 y ears, probably t he av erage for a 
classical nova (e.g., lOrio et all 120011) . In IShara et al.l 
([197S), Tcff of the constant bolometric luminosity phase 
is 780 000 K for a 1.25 Mq WD. Tes of the calcula- 
tions is not always published in pape r s on nova-outburst 
models (cf. IKovetz fc PrialnikI 119941: IPrialnik fc Kovetd 
1995), but published valu es range from 460 000 K for a 
0.8 M© (jShara et al.|[l979l ) to an extreme 2 500 000 K for 
a 1.4 M0 WD accreting at very high rate. The value of 
Tcff is dependent on mass-accretion rate, on the chemi- 
cal composition and on the WD temperature at the on- 
set of hydrogen burning, but generally a value of about 
700 000 K is only reached for M > 1.1 (Prialnik 
2009, priv. comm.). To summarize, the peak Tcff of 
about 740 kK indicates that the WD is very massive, with 
« 1.1 - 1.2 Mq. 

To explain the apparently decreased nitrogen abun- 
dance in the months after the outburst fFigure fTU]) . we 
conjecture that probably new material is already been 
accreted at this stage, while CNO burning at the bot- 
tom of the envelope is already proceeding at a lower rate 
before turn-off. 

A crucial point for the understanding of processes dur- 
ing the nova outburst may be the time-dependent predic- 
tion of surface abundances as well as the spectral analysis 
of high-resolution X-ray spectra taken from outburst to 
the end of surface H burning. The inspection of available 
Chandra and XMM-Newton observations of V4743 Sgr, 
taken in the 18 months after the outburst has shown that 
obtaining new X-ray grating spectra of future novae ev- 
ery week is highly desirable. 

In applications to novae, TMAP still lacks some 
physics, most notably the velocity fields. Although 
TMAP is not especially designed for the calculation of 
SEDs at extreme photospheric parameters, it is a flex- 
ible and robust tool for the determination of basic pa- 
rame ters like Tcff for line identifications (cf. IRauch et al.l 
|2008[) . and to derive a reliable range of abundances in 
the white-dwarf atmosphere. 
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